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Abstract 

Background:T\r\e recently solved solution structure of HCV (hepatitis C virus) p7 ion cliannel provides a solid structure basis 
for drug design against HCV infection. In tlie p7 channel the ligand amantadine (or rimantadine) was determined in a 
hydrophobic pocket. However the pharmocophore (— NH2) of the ligand was not assigned a specific binding site. 

Resu/ts:The possible binding sites for amino group of adamantane derivatives is studied based on the NMR structure of p7 
channel using QM calculation and molecular modeling. In the hydrophobic cavity and nearby three possible binding sites 
are proposed: His! 7, Phe20, and Trp21. The ligand binding energies at the three binding sites are studied using high level 
QM method CCSD(T)/6-31 1-I-G(d,p) and AutoDock calculations, and the interaction details are analyzed. The potential 
application of the binding sites for rational inhibitor design are discussed. 

Conclusions: Some useful viewpoints are concluded as follows. (1) The amino group (— NH2) of adamantane derivatives is 
protonated (— NHb^), and the positively charged cation may form cation-7t interactions with aromatic amino acids. (2) The 
aromatic amino acids (Hisl7, Phe20, and Trp21) are the possible binding sites for the protonated amino group (— NHb^) of 
adamantane derivatives, and the cation-7t bond energies are 3 to 5 times stronger than the energies of common hydrogen 
bonds. (3) The higher inhibition potent of rimantadine than amantadine probably because of its higher pKg value 
(pKa= 10.40) and the higher positive charge in the amino group. The potential application of p7 channel structure for 
inhibitor design is discussed. 
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Introduction 

Hepatitis C is an infectious disease [1] caused by the hepatitis C 
virus (HCV) [2], primarily affecting the liver. Hepatitis C is the 
leading cause for liver diseases in the USA. About 200 million 
people are infected with HCV worldwide [3]. The chronic 
infection of HCV can lead to scarring of the liver and ultimately to 
cirrhosis, which is generally apparent after many years. In some 
cases, those with cirrhosis will go on to develop liver failure, liver 
cancer or life-threatening esophageal and gastric varices [1]. So far 
there is no clinically proven vaccine [4,5] , and the most common 
therapy is based on a combination therapy of pegylated interferon- 
alpha (PEG-IFNa) and ribavirin (RBV), which only has a success 
rate of around 50% as well as severe side effects [6,7]. 
Development of more effective new drugs is absolutely necessary. 

The p7 channel plays multiple roles in virus life cycle and has 
several biological functions in HCV infection. Therefore the HCV 
p7 protein has been sought after as a potential anti-HCV drug 



target [8,9]. The p7 is a 63-residue membrane protein that 
oligomerizes to form ion channels with cation selectivity, for Na""", 
K"*", and Ca^"*" [10-13], and a more recent study has reported that 
the p7 channel mediated H intracellular conductance [14]. The 
adamantane derivatives and other several compounds [10—15] 
have been used in HCV clinical trials, but large variation in drug 
efficacy among the various HCV genotypes has been difficult to 
explain, because the drug target structure information of 
adamantane derivatives was not available, and the drug-target 
interaction mechanism was not clear. 

Recently the p7 channel structure in solution was solved by 
Chou and colleagues [16] using the state of the art NMR 
techniques. The top view and bottom view of the p7 channel are 
shown in Fig. 1 (A) and (B), respectively. In the p7 channel there 
are six equivalent hydrophobic pockets between the peripheral 
and pore-forming helices, consisting of Leu 52, Val 53 Leu55 and 
Leu 56 from H3, and Phe 20, Val 25 and Val 26 from H2 [16]. In 
Fig. 1 these hydrophobic residues are shown in green lines. The 



PLOS ONE I www.plosone.org 



1 



April 2014 I Volume 9 | Issue 4 | e93613 



Binding Sites in HCV p7 Channel 



ligand amantadine (or rimantadine) is located in the hydrophobic 
cavities. Fig. 2 is a close-up view of the binding location of ligand 
amantadine in the p7 ion channel, which is drawn based on the 
description in ref [16]. The binding location of amantadine (or 
rimantadine) in p7 ion channel, identified by Chou's lab, is 
different from the binding mode in the influenza M2 channel [1 7- 
19]. It turns out to be crucial to the functioning of the drug 
mechanism. While the tiny M2 channel in influenza gets plugged 
up by the drug molecule [18], in the p7 channel the drug nestles 
into a series of pockets within a folded outside edge of the funnel. 
When the drugs are in those pockets, the channel is unable to 
"exhale" and thus release ions [16]. 

In Chou's binding model of the p7 ion channel the amantadine 
(or rimantadine) is located in the hydrophobic cavity comprised by 
Phe20, val25, val26, Leu52, Val53, Leu55, and Leu56. The 
hydrophobic adamantane body of amantadine is surrounded by 
above hydrophobic amino acids in the cavity. However, the 
binding site of the adamantane's pharmocophore (— NH2) in the 
p7 ion channel was not accurately reported in ref [16], as saying 
the amino group (— NH2) of amantadine "on average points to the 
channel lumen" [16]. 

A pharmacophore is an abstract description of molecular 
features of a drug family, which are necessary for molecular 
recognition of a ligand by a biological macromolecule. Usually the 
drug pharmocophore makes the main contribution to the binding 
free energy and plays the key role in inhibition activity. The 




Figure 1. NMR solution structure (PDB code: 2iVI6X) of HCV p7 
ion channel and the positions of ligand amantadine. In p7 Ion 

cfiannel tfiere are six equivalent hiydrophobic pockets between tfie 
peripfieral and pore-forming helices. The ligand amantadine (or 
rimantadine) is located in the hydrophobic cavities. The pocket consists 
of Phe 20, Val 25, Val26, Leu52, Val53, Leu55, and Leu56,. The amino 
group of amantadine on average points to the channel lumen [1 6]. (A) 
A view from the top of channel. (B) A view from the bottom of channel. 
doi:1 0.1 371/journal.pone.009361 3.g001 




Figure 2. A close view of the binding pocket of amantadine in 
the p7 ion channel. The hydrophobic residues (Phe20, Val25, Val26, 
Leu52, Val53, Leu55, and Leu56) are shown in green line drawing, which 
comprise the hydrophobic binding pocket of amantadine. The positions 
of three possible binding sites Hisl7, Phe20, and Trp21 for the 
protonated pharmocophore group (-NH3*) of amantadine are shown 
in yellow, light blue, and white, respectively. All three aromatic residues 
(His17, Phe20, and Trp21) are on the chain 2. 
doi:10.1 371/journal.pone.009361 3.g002 

lUPAC defines a pharmacophore to be "an ensemble of stem and 
electronic features that is necessary to ensure the optimal supramolecular 
interactions with a specific biological target and to trigger (or block) its 
biological response" [39]. In the M2 proton channel [17-19] of 
influenza A virus the ligand adamantane derivatives may bind at 
more than one position, or the ligand binding location may change 
at deferent biological stages [20]. 

The authors of ref [16] pointed 'Hhe relatively poor stability of the 
protein— drug complex at the current stage of our study precludes fill-scale 
structure determination" . In this situation molecule modeling and 
accurate QM calculation may help to solve the binding sites of the 
pharmocophore group (— NH2) within the hydrophobic pocket of 
p7 ion channel. 

Method and Theory 

Adamantane derivatives are alkaline compounds. The pK^ 
values of residues are important for binding interactions. In 
proteins the pK^ values of residues can change in a broad range 
because of the influence from the interaction environment. The 
pKj, value of an ionizable amino acid is evaluated using the 
following equation, 

pK„=pK^°^ + ApK:r (1) 

In Eq. 1 the pK;,™"'^ is a model value assigned to certain amino 
acid types. For histidine the model value pKj;™"'' = 6.50. The term 
ApKa'^"^ is the correction value from the environment residues in 
protein. The acidic dissociation constants pK^i values of all acidic 
and alkaline residues in p7 ion channel are calculated using 
PROPKA3.1 software package [21-24]. 

The binding free energies between ligand and host protein at 
three possible binding sites (His 17, Phe20, and Trp21) are 
calculated by using AutoDock 4.0 software package [25], and 
the grid maps of the protein used for docking process were 
calculated with the AutoGrid 4.0 [26]. The grid dimensions 
centered at the defined pocket were 40x40x40 grid points with a 
spacing of 0.375 A in each dimension. Gasteiger charges, 
computed by ADT (AutoDock tools), were assigned to both 
amantadine and receptor. Docking simulations were performed 
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with the Lamarckian Genetic Algorithm [25-27], using maximum 
number of 250,000 energy evaluations, mutation rate of 0.02, 
cross over rate of 0.08, and elitism value of 1 . All other docking 
parameters were left at the default values. Each docking job 
included 200 independent runs. Finally, the docked poses for each 
ligand within 1.0 A in the root mean square deviation (rmsd) 
tolerance of each other were clustered together with the really 
close binding energy. In the calculations a flexible model is used 
for both ligand and host acceptor. It means that the Kgand and 
host protein can adjust their conformations and orientation to 
make the best docking effect. 

The protonated amino group (— NHg""") of amantadine (or 
rimantadine) can forms a stable cation-7t bond with the aromatic 
residues. The physical nature and properties of cation-7t interac- 
tions [28-32] are very dilferent from the well known three type 
interactions: van der Waals interaction, electrostatic interaction, 
and hydrogen bond interaction. Some available force field 
parameters may not include the cation-7t interactions correctly. 
In cation-7i interactions the electron correlation and dispersion 
interaction make important contribution, which have to be 
described by using post Hartree-Fock methods, such as MP2, 
CCDS or CCSD(T) [33-36]. In this study the structures of 
molecular monomers are optimized using CCSD/6-31 1+G(d,p), 
and the cation-i: interaction energies are calculated at CCSD(T)/ 
6-311+G(d,p) level. 




Figure 3. The ligand-acceptor binding interaction between 
amantadine and p7 ion channel at binding site His17. (A) The 

docking structure between amantadine and p7 channel at the site 
His17. (B). QM calculation for cation-7t interaction between CH3NH3* 
and the aromatic side chain of His. The protonated amino group (— 
NH3*) of amantadine perpendicularly points to the 71-plane of His17, 
and forms a stable cation-71 bond. The cation-jt bond length is 3.072 A 
and the interaction energy is -50.28 kJ/mol. 
doi:1 0.1 371/journal.pone.009361 3.g003 



Results 

In this section the computational and modeling results are 
reported using tables and figures, and followed by brief illustration 
and comparison. 

Acidic Dissociation Constants (pKa) of Residues 

The pharmocophore (— NHj) of rimantadine and amantadine 
are alkaline group that is sensitive to the acidic dissociation 
constants (pKJ values of amino acids. The pK^ values of residues 
in proteins can change in a broad range because of the influence 
from the interaction environment. The pK^ values of residues in 
p7 ion channel are calculated by using the PROPKA3.1 software 
[22-24,37], and the results are listed in Table 1. 

Among the 378 amino acids (63x6) only 18 histidine residues 
(His 17, His37, His59) are acidic, and the pK^i values of histidine 
residues are in the region 5.44 to 6.22. On the other hand, 42 
residues (Lys3, Lys33, Trp42, Arg35, Arg54, Arg57, and Arg60) 
are alkaline amino acids. Among them the pKj, values of 24 
arginine residues are larger than pK3>12. The p7 cation channel 
is an alkaline amino acid dominated protein. 

In the last two lines of Table 1 the pK^i values of amantadine 
(pKa = 9.00) and rimantadine (pR, = 10.40) are checked from ref 
[38]. In common cell condition most molecules of amantadine and 
rimantadine appear in protonated form (Ad— NH^"^. 

Possible Binding Sites for Amantadine Pharmocophore 

(-NH^) 

As described in ref [16], the binding location of amantadine (or 
rimantadine) in p7 ion channel is in the hydrophobic pocket 




Figure 4. The ligand-acceptor binding interaction between 
amantadine and p7 channel at binding site Phe20. (A) The 

docking structure between amantadine and p7 channel at the site 
Phe20. (B). QIVI calculation for cation-71 interaction between CH3NH3* 
and the aromatic side chain of Phe. The protonated amino group (— 
NH3*) of amantadine perpendicularly points to the 71-plane of Phe20, 
and forms a stable cation-71 bond. The cation-7t bond length is 3.084 A 
and the interaction energy is -62.80 kJ/mol. 
doi:1 0.1 371/journal.pone.009361 3.g004 
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Table 2. Binding free energies and catlon-ji Interaction energies of annantadlne (or rimantidine) at three possible binding sites 
(Hisl7, Phe20, and Trp21) of p7 Ion channel. 





*Docking calculation: 


Binding sites 


Hisl7 


Phe20 


Trp21 


Free energy 


-19.09 


-22.81 


-24.53 


**CCSD(T)/6-31 1+G(d,p) calculation: 


Interaction pair 


CH3NH3+-His 


CHsNHs^-Piie 


CH3NH3+-Trp 


Cation-Ti energy 


-50.278 


-62.80 


-82.53 



*Calculated using Autodock4.0 [33], energy in kJ/mol. 

**CCSD/6-311G{d,p) [37-40], energy in l<J/mol. In QM calculations the protonated amantadine is simplified as CH3NH3^, and the aromatic amino acids are replace by 

their aromatic side chains. 

doi:l 0.1 371 /journal.pone.009361 3.t002 



comprised by Phe20, val25, val26, Leu52, Val53, Leu55, and 
Leu56. This result is overall consistent with a mutational study 
[38] showing that mutations in residues 50-55 significantly reduce 
drug sensitivity of the channel. 

The pharmocophore (— NH''"3) of amantadine is an alkaline 
group, also a cation. The possible binding interactions are acid- 
alkaline interaction and cation-7t interaction. In the hydrophobic 
pocket and nearby there are three aromatic amino acids (His 17, 
Phe20, and Trp21). The His 17 is also an acidic residue. The three 
residues are the possible binding sites for the protonated amino 
group (— NHs""^ of amantadine and rimantadine. The positions of 




Figure 5. The ligand-acceptor binding interaction between 
amantadine and p7 channel at binding site Trp21. (A) The 

docking structure between amantadine and p7 channel at the site 
Trp21. (B). QM calculation for catlon-n interaction between CH3NH3* 
and the aromatic side chain of Trp. The protonated amino group (— 
NH3*) of amantadine perpendicularly points to the n-plane of Trp21, 
and forms a stable catlon-jt bond. The cation-71 bond length Is 2.99 A 
and the Interaction energy Is -82.53 kJ/mol. 
dol:1 0.1 371/journal.pone.009361 3.g005 



three residues (Hisl7, Phe20, and Trp21) are shown in Fig. 2. All 
the three residues are on the chain 2. 

(1) Binding site Hisl7. Histidine is the only acidic amino 
acid in p7 protein, possessing the pK^j = 5.44~6.22 in p7 channel, 
littie lower than 7. In the common cell condition histidine may 
appear in either protonated form or neutral form. When the 
residue His 17 is in neural form, the protonated amino group (— 
NHj""^ of amantadine (or rimantadine) can forms a stable cation-71 
bond [28-32] with the aromatic side chain imidazole of Hisl7. A 
docking structure of cation-7i interaction between amantadine and 
His 17 is shown in Fig. 3 (A), and Fig. 3 (B) is an illustration 
portrait for the cation-71 interaction between CH3NH3''" and 
imidazole (side chain of His 17), in which the protonated amino 
group (— NHj"^ perpendicularly point to the 7:-plane of imidazole 
ring. The distance from nitrogen atom to the Tt-plane is 3.1 A, and 
the bond energy is —50.28 kJ/mol, 2.5 times of the common 
hydrogen bond energies ( 20 kJ/mol). 

Among the 20 natural amino acids histidine is the unique 
member that can join 5 type molecular interactions (hydrogen 
bond, cation-7t, polar hydrogen-7i, 7i-7t stocking, and coordinate 
interaction) [28] . The pK^j value of histidine can change in a broad 
range around 6. 0, affected by the interaction condition. 
Consequendy histidine could be a proton donor or acceptor, 
and appear in neutral form or protonated form. In the p7 ion 
channel the His 17 residues may play an important role in the 
cation conductance. Therefore the binding of amantadine on the 
His 17 may significantly affect the biological function of p7 ion 
channel. 

(2) Binding site Phe20. In the p7 ion channel [16] the Phe20 
is an artificial mutation that replaced the Leu20 in wUd HCV p7 
channel. The aromatic amino acid Phe20 is a part of the 
hydrophobic pocket and directly contacts with the hgand 
amantadine or rimantadine [16], as shown in Fig. 2 (A). The 
distance between Phe20 and amino group of amantadine is 3.6 to 
4.0 A, in the cation-7i interaction region. Fig. 4 (A) shows the 
docking structure between amantadine and Plie20, in which the 
protonated amino group (— NHs"^ perpendicularly points to the 7i- 
plane of Phe20, and a stable cation-7t bond is formed. Fig. 4 (B) is 
an illustration portrait for the cation-7t interaction between 
CHaNHs^ and CgHg (side chain of Phe20). The optimized 
distance from nitrogen atom to the 7t-plane is 3. 1 A, and the 
cation-71 interaction energy is —62.80 kJ/mol, three folds as the 
common hydrogen bond energy. 

(3) Binding site Trp21. The tryptophan is the largest 
aromatic amino acid in the 20 natural amino acids. In the p7 
channel the positions of Trp2 1 and Phe20 are on the two sides of 
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the ligand amantadine: Trp21 is close the amino group, and 
Phe20 is near the adamantine body of the amantadine. The side 
chain of Trp2 1 plugs into the p7 channel. Similar to the Phe20, 
the aromatic amino acid Trp2 1 can form cation-7t interaction with 
the protonated amino group (— NH^^) of amantadine (or 
rimantadine). The binding interaction structure of amantadine- 
Trp2 1 is shown in Fig. 5 (A). The cation-K interaction energy (— 
82.53 kj/mol) between Trp21 and amantadine is four folds as the 
common hydrogen bond energy, which is the largest one among 
the three aromatic amino acids because of its large aromatic 7t- 
system. 

Calculations of Binding Free Energies 

In the above calculations the binding energies between 
amantadine and three aromatic amino acids are calculated by 
using QIV'I method. The binding energies are the main contribu- 
tion to the binding free energies. However, the full binding free 
energies have to be calculated using docking method. The docking 
free energies between ligands amantadine and p7 channel at the 
three binding sites (Hisl7, Phe20, and Trp21) are calculated using 
AutoDock4.0 [25]. In the calculations a flexible model [26,27] for 
both ligand and host acceptor is used, meaning that the local 
conformations of both ligand and host protein can be adjusted to 
reach the best docking effect. The calculated docking free energies 
are listed in Table 2. 

In Table 2 the binding free energies of docking calculations are 
around —20 kJ/mol, much smaller than the cation-7t interaction 
energies that are [:alculated using high level quantum chemical 
method CCSD(T)/6-311+G(d,p) [33-36]. The binding free 
energies include all interaction terms between ligand and host 
protein. The cation-7t interaction energy is the main contribution 
to the binding free energy. In the docking calculations the binding 
free energies depend on the force field parameters [30,40-43]. 
However, the cation-7i interaction energies may be not correcdy 
described by the available force field parameters [30,40-43]. 
Despite the value differences between binding free energies (dock 
CEilculations) and cation-7i energies (QM calculations), the order of 
the two calculation methods are the same: Trp21>Phe20>HisI7. 
As shown in Table 2, the energies of cation-7t interactions are 
much larger (2.5 to 4 folds) than the energies (- — 20 kJ/mol) of 
common hydrogen bonds. 

Discussion 

In this study the binding-site searching for pharmocophore 
group (— NH^) of adamantane derivatives are performed in the 
hydrophobic binding pocket comprised by Phe20, val25, val26, 
Leu52, Val53, Leu55, and Leu56, which is determined by NMR 
experiment [16], and supported by the previous mutation 
experiments [15,38,46]. It is common that the ligands may have 
more than one biding sites [44] in the host proteins, and the ligand 
binding sites may change at different biological stages [45] . Within 
the binding pocket of adamantane derivatives, three possible 
binding sites (His 17, Phe20, and Trp21) for the amino group of 
amantadine (or rimantadine) are suggested based on high level 
QM calculations and docking calculations. 

In the binding energy calculations we pay high attention on the 
acid-alkaline interaction and cation-7t interaction, because the 
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Conclusion 

In the paper [16] by Chou and colleagues, little attention was 
given to the pharmacophore amino group of amantadine and 
rimantadine, so this is now considered in this current manuscript. 
This study focuses on the recently published NMR structure of 
HCV p7 and aims to analyze this structure with respect it's 
potential binding to adamantine derivatives using molecular 
modeling and quantum mechanics. Based on the modeling results 
three potential binding sites (His 17, Phe20 and Trp21) are 
proposed. The best binding site is the Trp21. However, the 
protonated amino group might change the original binding site to 
Phe20, which could explain the resistance-conferring Leu20Phe 
mutation. 

All of these residues have been proposed before and have been 
characterized by mutagenesis, so the data are in line with results 
from experiments performed in biological systems. Insights are 
also provided with respect to the know Leu20Phe resistance 
mutation. The results regarding the binding location of amanta- 
dine or rimantadine are confirmatory of the results described in 
the paper [16] by Chou and colleagues. However, in the current 
manuscript further results are presented regarding the potential 
binding of the protonated amino group. We hope that the possible 
binding sites and the binding interactions, proposed in this study, 
may help the HCV inhibitor design targeting the p7 channel 
protein. 
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